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Abstract. We review our theoretical approach to neutral current photon emission on nucleons and nuclei in the few-GeV 
energy region, relevant for neutrino oscillation experiments. These reactions are dominated hy the weak excitation of the 
A(1232) resonance but there are also important non-resonant contributions. We have also included terms mediated by nucleon 
excitations from the second resonance region. On nuclei, Pauli blocking, Fermi motion and the in-medium A resonance 
broadening have been taken into account for both incoherent and coherent reaction channels. With this model, the number 
and distributions of photon events at the MiniBooNE and T2K experiments have been obtained. We have also compared to 
the NOMAD upper limit at higher energies. The implications of our findings and future perspectives are discussed. 
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INTRODUCTION 

A good understanding of (anti)neutrino cross sections is crucial to reduce systematic uncertainties in oscillation 
experiments [1], Our present knowledge of neutrino-nucleus interactions has been improved by a new generation 
of oscillation and cross section experiments. Over the last decade, K2K, NOMAD, MiniBooNE, SciBooNE, MINOS, 
and more recently T2K and MINERvA have obtained a wealth of data on quasielastic(-like) scattering, incoherent and 
coherent single pion production, and inclusive cross sections. These results challenge our understanding of neutrino 
interactions with matter and have triggered a renewed theoretical interest [ 2 ]. 

One of the possible reaction channels is photon emission induced by neutral cutTent (NC) interactions (NC 7 ), which 
can occur on single nucleons and on nuclear targets with incoherent or coherent reaction mechanisms. Weak photon 
emission has a small cross section compared, for example, with pion production, the most important inelastic reaction 
channel. In spite of this, NC photon emission turns out to be one of the largest backgrounds in —?> Vg (v^ ^ w) 
oscillation experiments where electromagnetic showers from electrons (positrons) and photons are not distinguishable. 

The first effort to put the description of NC photon emission on solid theoretical grounds was reported in Ref. [3]. 
The reaction on nucleons was studied with a microscopic model developed in terms of hadronic degrees of freedom: 
nucleon, A(1232) resonance and mesons. Coherent photon emission off nuclear targets was also evaluated, treating 
the nucleus as a scalar particle and introducing a form factor to ensure that the coherence is restricted to low- 
momentum transfers. The NC 7 reactions on nucleons and nuclei were also studied using a chiral effective field theory 
of nuclei [4, 5, 6 ], phenomenologically extended to the intermediate energies (Ey 1 GeV) in Ref. [7]. In this model, 
a rather strong in-medium suppression of the A(1232) excitation is compensated by rapidly growing contact terms. 

Our theoretical model [ 8 ], summarized below, extends and improves several relevant aspects of the existing descrip¬ 
tions. The energy and target mass dependence of the integrated cross sections are displayed and discussed. The model 
is then applied to the specific conditions of the MiniBooNE, T2K and NOMAD experiments. 


FORMALISM 
NCy on nucleons 


The cross section for 


V, v{k) + N{p) v,v{k') + N{p') + Y[ky) 


( 1 ) 


in the laboratory frame is given by 


d^<^{V,V) _ Ey r d^k' ly^v) ng 

dEydD.{ky) ~ \k\ 1671^ J \k'\ 

in terms of contraction of the leptonic tensor 
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with the hadronic one 
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In these equations, M denotes the nucleon mass, Ey and E'^ are the photon and final nucleon laboratory energies 
and q — k — k' is the 4-momentum transfer. The model is defined by the set of Feynman diagrams for the hadronic 
current [8]. 

{Nrlj^cyi^m = a{p')r^puip)e;iky), (5) 

shown in Fig. 1; e{ky) is the photon polarization vector. 



FIGURE 1. Feynman diagrams for the hadronic current of NC photon emission considered in Ref. [8]. The first two diagrams 
stand for direct and crossed baryon pole terms with nucleons and resonances in the intermediate state: BP and CBP with B = N, 
A(1232), A(1440), A(1520), A(1535). The third diagram represents the r-channel pion exchange: TtEx. 

The amputated amplitude F^^ for nucleon pole terms, NP and CNP, is given by 
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with currents 


4c(^) = fFii‘l^) + :^<y^^qpF2{q^)-fr5FA{q^), (7) 

■fEMih) = rF,{0) + :^a^^iky)yF2{0). ( 8 ) 


At threshold, this mechanism is fully constrained by gauge and chiral symmetries, and the partial conservation of 
the axial current (PCAC). These terms are infrared divergent when the photon energy Ey ^ 0 but this doesn’t have 
practical consequences because very soft photons will not be experimentally resolved. The extension towards higher 
energy transfers required to make predictions at Fy ~ 1 GeV is performed using phenomenological parametrizations 
of the weak form factors. Isospin symmetry allows to relate the NC vector form factors 


f/J = (l-4sin2 0 ^)Fg-FW-FW 

Fg = (l-4sin2 0 ^)Fg-Fg-Fg, 


to the electromagnetic ones that have been extracted from electron scattering data. For the axial form factors 
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we have adopted a conventional dipole parametrization [9]. Strange form factors, whose present values are consistent 
with zero, have been neglected. 

The AP and CAP terms 
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with nucleon-A(1232) transition currents 
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is given in terms of vector 

Ctsiq^) = (1 -2sin2 Bw)Cl_^{q^) (15) 

and axial {C^_^ = ( 2 ^ 5 ) form factors. The vector form factors are related to the helicity amplitudes extracted in the 
analysis of pion photo- and electro-production data. We have adopted the parametrizations of the helicity amplitudes 
obtained with the unitary isobar model MAID [10]. Assuming the so called Adler model [11, 12], the axial current 
depends on just one form factor 

Cj(,^)=c?(0)(l-j|-) 7 (16) 

For it we have adopted a dipole ansatz with Cf(0) = 1.00 ±0.11 and Ma = 0.93 GeV, following the fit to v^c/ 
p-A++n BNL and ANL data [13]. 

We have extended the model to higher energy transfers by taking into account intermediate states from the second 
resonance region, which includes three isospin 1/2 baryon resonances Pii(1440), Di3(1520) and 5ii(1535). The 
structure of the contribution of Pii(1440) and 5ii(1535) to the amputated amplitudes is similar to the one of the 
nucleon, while the A(1520) amplitudes resemble the A(1232) ones. Detailed expressions can be found in Ref. [ 8 ]. 
With these N* terms we have adopted a strategy similar to the one for the A(1232) terms. The vector form factors are 
expressed in terms of the empirical helicity amplitudes extracted in the MAID analysis. As there is no experimental 
information to constrain the axial form factors, following Ref. [14], we have kept only the leading axial terms and 
used PCAC to derive off-diagonal Goldberger-Treiman relations between the corresponding axial couplings and the 
N* -A- Nn partial decay widths. For the dependence we have assumed a dipole ansatz like in Eq. (16) with a natural 
value of M^ = 1.0 GeV. 

Finally, the nEx mechanism originates from the Z 7 ;r*' vertex fixed by the axial anomaly of QCD 
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It is nominally of higher order [4] and, indeed, gives a very small contribution to the cross section as will be shown 
below. There are other f-channel p and (O exchange mechanisms [3] that arise from the anomaly-mediated Z'^yp and 
Z^^yo) interactions [15]. Among them, the O) contribution is favored by the size of the couplings. Reference [4] assumes 
that the p and (O exchange mechanisms, taken from Ref. [3], saturate the low-energy constants in the contact terms. 
The contribution of these contact terms to the NC 7 cross section on the nucleon is very small at Ey < 550 MeV [3, 4], 
as expected from power counting arguments. The extension to higher energies requires the introduction of poorly 
understood form factors [3, 7]. The cross section from these mechanisms increases fast with energy, which might be a 
concern for experiments at high energies, or with a high-energy tail in the neutrino flux (like T2K). Nevertheless one 
should recall that this trend will be limited by unitarity bounds; these amplitudes will be modified by loop contributions 
and partially canceled by contact terms of even higher orders. 



NCy on nuclei 


This reaction on nuclear targets can be incoherent 


V, v(^) + Az-^ v,v{k') + Y[ky) + X 


( 18 ) 


or coherent 


v,v{k) + Az\gs{pA) v,v{k') +Az\g,{p'A) + 7{ky) 


(19) 


depending on whether the final nucleus is in an excited or in its ground state. 

In a many-body scheme adapted to (semi)inclusive reactions on finite nuclei by means of the local density approxi¬ 
mation, the NC 7 incoherent cross section is 



( 20 ) 


The hadronic tensor is obtained from the imaginary part of the contributions to the selfenergy with a single photon 
in the intermediate state. In a density expansion, the lowest order contribution is depicted in Fig. 2. 
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FIGURE 2. Diagrammatic representation of the one-particle-one-hole-photon (Iplhy) contributions to the self-energy in 
nuclear matter. The black dots represent Z^N —>■ yN elementary amplitudes. To obtain the imaginary part, the intermediate states 
intersected by the dashed line have to be placed on the mass shell. 


After the simplification of evaluating the F^^ amplitudes at an average nucleon hole four momentum (p^), 
becomes (see Ref. [ 8 ] and references therein for details) 



( 21 ) 


where 



( 22 ) 


while UR{q — ky,kp,kp) stands for the Lindhard function (definition and explicit expressions can be found in 


Ref. [16]). The Fermi momentum depends on the local density of nucleons in the nucleus via kp{r) = [ 37 r^pAr(r)]*/^. 
The density distributions are based on empirical determinations in the case of protons and on realistic theoretical 


models in the case of neutrons. 

It is known that the A(1232) properties are strongly modified in the nuclear medium. This important nuclear 
correction is taken into account by the following substitutions in the A(1232) propagator 


Ma —>■ Mjg + Rel^A (p), 
Ta ->• rA-2Im5:A(p). 


(23) 

(24) 


The real part of the in-medium A selfenergy, Ea, receives an attractive (negative) contribution from the nuclear mean 
field, which is partially cancelled by an effective repulsive piece from iterated A-hole excitations. As the net effect 
is smaller than the precision achievable in current neutrino experiment we simply take ReEA(p) ~ 0. The resonance 
decay width is reduced to Fa because the final nucleon in A nN can be Pauli blocked but, on the other hand, it 
increases because of the presence of many body processes such as AN NN, AN NN n and ANN NNN. 
These new decay channels, accounted in ImEA, have been parametrized as a function of the local density in Ref. [17]. 

For the coherent reaction of Eq. (19) one has [ 8 , 18] that 
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where stand for the amputated photon production amplitudes for the different mechanisms i = 

NP,CNP, nEx,RP,CRP [R — A,N{lA-AQ),N{\5?i5),N{\520)]. In this case, the nucleon wave functions remain 
unchanged so that one has to sum over all nucleons at the amplitude (not amplitude squared) level. This leads to 
the trace in Eq. (28) and the nuclear density distributions in Eq. (27). Therefore, the coherent production process is 
sensitive to the Eourier transform of the nuclear density. In the elementary Z^N Ny process, energy conservation 
is accomplished by imposing q^^ = Ey, which is justified by the large nucleus mass. The transferred momentum is 
assumed to be equally shared between the initial and final nucleons. A lengthy discussion about this prescription and 
the local treatment of the A propagation can be found in Ref. [8]. The modification of the A in the medium outlined 
above has also been considered here. 


NCy CROSS SECTIONS 
On nucleons 

The integrated NCy cross sections on protons and neutrons as a function of the (anti)neutrino laboratory energy are 
displayed in Eig. 3. As in other processes, the different helicities of v and v are responsible for different interferences, 
resulting in smaller v cross sections with a more linear energy dependence. The A mechanism is dominant. Its 
contribution is the same on protons and neutrons because of the isovector nature of the electroweak A — A transition. 
The contribution of NP + CNP terms is also important, being only about 2.5 smaller than the A one at Ey)!/) ~ 1 -5 
GeV. Above ~ 1.5 GeV, the A(1520) contribution is also sizable and, for v on protons, comparable to the one of 
NP + CNP. However, the A(1440) and A(1535) and nEx mechanisms can be safely neglected. The fact that the 
A(1520) resonance is the only one, besides the A(1232), playing a significant role at Ey < 2 GeV has also been 
observed in pion production [19] and for the inclusive cross section [14]. More details, including differential photon 
energy and angular distributions as well as comparisons with earlier studies can be found in Ref. [8]. 


On nuclei 

In the left panel of Eig. 4, we show our predictions for the (anti)neutrino incoherent photon emission cross sections 
on as a function of the (anti)neutrino energy. It is clear that neglecting nuclear medium corrections is a poor 
approximation. By taking into account Eermi motion and Pauli blocking, the cross section already goes down by more 
than 10%. With the full model that also includes the A resonance in-medium modification, the reduction is of the order 
of 30%. The right panel of Eig. 4 shows total NCy incoherent cross sections for different nuclei. The curves indicate 
an approximated A-scaling. Nevertheless, the cross section is smaller for heavier nuclei, particularly We should 

stress that the observed deviation from scaling cannot be explained only by neutron cross sections being smaller than 
proton ones. Photon distributions do not show any appreciable A dependence in the shapes as can be seen in Pig. 10 of 
Ref. [8]. Comparisons with other calculations are also displayed and discussed in that article. 

Integrated cross sections for coherent NCy are shown in Pig. 5. These are about a factor 10-15 smaller than the 
incoherent ones. Thus, their relative relevance is similar if not greater than in pion production. As is apparent from the 
left panel, the dominance of the AP + CAP terms is more pronounced in this case, with small corrections from A( 1520) 
excitation. Nucleon-pole contributions are negligible because the coherent kinematics favors a strong cancellation 
between the direct and crossed terms. The nEx terms vanish exactly for isospin symmetric nuclei because amplitudes 
for protons and neutrons cancel with each other. 
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FIGURE 3. vN —>• vNy (left panel) and vN vNj (right panel) cross sections on protons and neutrons [8]. A cut of Ey > 
140 MeV in the phase space integrals has been applied. The error bands in the full-model results (solid lines) represent the 
uncertainty in the axial AA coupling C^(0) = 1.00 ±0.11 according to the determination of Ref. [13]. The curves labeled N, 
A, Di 3 , Pii and An stand for the partial contributions of the BP and CBP mechanisms of Fig 1; label Tt corresponds to kEx. The 
N* contributions have been removed from the results labeled “no A*”. 


Unlike it and p f-exchange terms, the O) contribution does not vanish for symmetric nuclei because amplitudes on 
protons and neutrons add up. In Ref. [3] it was found that the coherent NC 7 from co exchange plays a sub-dominant 
role at fiv 1 GeV, compared to naive estimates, being suppressed by form factors and recoil. On the other hand, 
because of their strong energy dependence of the contact terms, the coherent cross section in Ref. [7] is dominated by 
contact terms for Ey y > 650 MeV. However, as discussed after Eq. (17), these results are not only highly sensitive to 
unknown form factors but should also be constrained by unitarity. 

Finally, the right panel reveals that the nuclear dependence is stronger than for the incoherent channel. The coherent 
cross sections neither scale with A, like the incoherent one approximately does, nor with as one would expect from 
the dominant isoscalar AP mechanism. This is related to the structure of the axial and vector currents atq^ ky, favored 
by the nuclear form factor: a more elaborated explanation can be found in Ref. [ 8 ]. 

SINGLE PHOTON EVENTS AT MINIBOONE, T2K AND NOMAD 
NCy and the low energy excess at MiniBooNE 

The MiniBooNE experiment, designed to explore the short-baseline Vg oscillations, has found an excess 

of electron-like events over the predicted background in both v and v modes [20, 21]. The excess is concentrated 
at 200 < Ey ^ < 475 MeV, where is the neutrino energy reconstructed assuming a charged-current quasielastic 
(CCQE) nature of the events. Recent studies indicate that explanation of this anomaly cannot reside in oscillations, 
even involving sterile neutrinos [22, 23]. On the other hand, it could have its origin in poorly understood backgrounds. 

At low the background is dominated by photon emission because Cherenkov detectors like MiniBooNE cannot 
distinguish electrons from single photons. The largest source of single photons is NCtt** production where one of 
the photons from the it^ 77 decay is not identified. This background has been constrained by the MiniBooNE’s 
NCtt*' measurement [24]. The second most important process is NC 7 . The MiniBooNE analysis estimated this 
background using the NCtt*' measurement, assuming that NC 7 events come from the radiative decay of weakly 

























FIGURE 4. Left panel: Incoherent NCycross sections on according to the model of Ref. [8]. All curves have been obtained 
with an Ey > 140 MeV cut in the phase space. Solid lines are from the complete model at the nucleon level, while the dotted 
lines were obtained without the N* contributions. Curves denoted as “Free” (upper blue curves) do not include any nuclear 
correction: (<Ja = ZOp +N<J„). Curves labeled as “Full” (lower red curves) take into account nuclear effects. The error bands 
show the uncertainty from the the axial AA coupling (C^(0) = 1.00it0.il). Right panel: Integrated cross sections for different 
nuclei (*^C,*®0,'*°Ar, ^®Ca,^®Fe and ^°*Pb) divided by the number of nucleons. 


produced resonances, mainly A —Ny [20, 21]. This procedure neither takes into account the existence of non-resonant 
terms, nor the coherent channel. If the NCy emission estimate were not sufficiently accurate, this would be relevant to 
track the origin of the observed excess. 

We have applied the model described above to calculate the number and distributions of single photon events 
at MiniBooNE [25], using the available information about the detector mass (806 tons) and composition (CH 2 ), 
the total number of protons on target (POT), 6.46 x 10^° in v mode and 11.27 x 10^® in v mode, [21], the flux 
prediction [26] and photon detection efficiency [27]. The yields from the incoherent channel are the largest ones. 
Those from the coherent channel and the reaction on protons, which are comparable, are smaller but signihcant. The 
coherent contribution is particularly important for antineutrinos and in the forward direction (see the plots in Ref. [25]). 

Our results for the distributions are shown in Fig. 6 . The error bands correspond to the uncertainty in C 5 (0). A 
more complete error analysis, leading to similar bands, can be found in Ref. [25]. The comparison with the MiniBooNE 
in situ estimate [21, 27] shows a good agreement: the shapes are similar and the peak positions coincide. The largest 
discrepancy is observed in the lowest energy bin. In the two bins with the largest number of events, the two calculations 
are consistent within our eiTorbars. For higher Ey ^ values, our results are systematically above the MiniBooNE 
estimate although the differences are very small. The error in the detection efficiency (~ 15%) [27], not considered in 
this comparison, will partially account for the discrepancies. The overall agreement is also good in comparison to the 
estimate of Zhang and Serot [7] in spite of the differences in the approaches, in contrast to the hndings of Hill [28], 
obtained with a rather high and energy independent detection efficiency, and neglecting nuclear effects. Therefore, we 
have found that neutral current photon emission from single-nucleon currents is insufficient to explain the event excess 
observed by MiniBooNE in both neutrino and antineutrino modes. 
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FIGURE 5. Left panel: Left panel: Coherent NCy cross sections on according to the model of Ref. [ 8 ]. A photon energy 
cut of Ey > 140 MeV has been applied. Red solid lines stand for results from the complete model, with error bands determined by 
the uncertainty of ±0.11 in C 5 (0) [13]. The solid blue lines below, labeled as “no N*’\ display the cross sections without the N* 
amplitudes, while the magenta dotted ones denote the contributions from the A mechanisms alone. Right panel: Integrated cross 
sections for different nuclei ('^C,*® 0 ,^*'Ar, ^®Ca,^®Fe and ^®^Pb) divided by the number of nucleons. 
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FIGURE 6. Ey® distributions of total NCy events for the v (left) and v (right) modes. The error bands correspond to the 
uncertainty in C 5 (0) = l.OOiO.ll [13] around the central value labeled as “1.0”. The curves labeled as “no N*” show results 
without the N* contributions. The “MB” histograms display the MiniBooNE estimates [27]. 


NCy events at the Super-Kamiokande detector 

The T2K experiment has measured Vg appearance in a beam, obtaining the first indication of a nonzero value 
of 013 [29]. The increasing precision allows for more detailed studies of neutrino properties and might lead to the 
discovery of CP violation in the lepton sector. Progress in this direction requires a better control over irreducible 
backgrounds. As MiniBooNE, the Super-Kamiokande (SK) water detector is a Cherenkov one and cannot distinguish 
photons from electrons. The largest part of this background comes from NCtt**, which has been significantly reduced 
using specific reconstruction algorithms [30], making the irreducible NC 7 contribution relatively more important. 























































With our model we have predicted the number and distributions of NCy events at the SK detector. We have used the 
flux of the off-axis neutrino beam from Tokai at SK [31] but neglect the tail above Ey = 3 GeV. For the recent T2K Ve 
appearance analysis, corresponding to A^pox = 6.57 x in v mode [30] we have obtained a total number of 

= 0.421 ±0.051, (29) 

before efficiency corrections, with the error corresponding to the uncertainty in C5(0). When compared to the 
equivalent calculation performed with the NEUT event generator [32], our result turns out to be 2-3 times larger. 
The disagreement is mostly in the normalization because the shapes of the photon energy and angular distributions 
are similar [33]. This finding is in line with the comparison of the NCy integrated cross sections on from different 
models presented in Fig. 9 of Ref. [34], where the NEUT result is below the rest. 


The NCy limit at NOMAD 


The high resolution data taken by the NOMAD experiment allow a sensitive search for neutrino induced single 
photon events. The experiment has obtained an upper limit of 4.0 x 10^^ single photon events per charged-current 
ones with 90 % CL, at Ey ~ 25 GeV [35]. Althogh the NCy models developed so far are not applicable at the high 
energy transfers that can occur in NOMAD, in the limited region of phase space where these models are valid, they 
should fulfil the NOMAD constraint as a necessary condition. In our case, restricting the invariant mass of the outgoing 
nucleon-photon pair to W < 1.6 GeV, where the model can be retained applicable, and neglecting nuclear effects (that 
would reduce the cross section) for simplicity, we obtain 


(7(NC7,W< 1.6GeV) 
(7(V^A fi^X) 


(30) 


at Ey = 25 GeV using NOMAD inclusive charged current cross section data [36] for the denominator. Our prediction 
is then safely below the NOMAD limit. 


OUTLOOK 

With the microscopic model developed in Ref. [8] we have calculated photon emission on nucleons and nuclei 
in a kinematic region of interest for current and future neutrino experiments. After the inclusion of N* excitation 
mechanisms from the three lightest states, the model can be considered reliable up to nucleon-photon invariant masses 
of around 1.6 GeV. The main uncertainty arises from the uncertainty in the dominant A —A(1232) axial coupling that 
can be constrained in weak pion production experiments. In the case of nuclear targets, we have also studied the mass 
dependence of both incoherent and coherent reaction channels emphazising the large (~ 30%) nuclear corrections. 

We have predicted single photon events at the MiniBooNE and SK detectors where they represent an important 
electron-like irreducible background component. In the case of MiniBooNE, our results are consistent with the in-situ 
estimate, obtained with a much poorer model. Based on this, we conclude that photon emission processes from single¬ 
nucleon currents cannot explain the excess of the signal-like events observed at MiniBooNE. As mentioned above. 



FIGURE 7. Diagrammatic representation of two-particle-two-hole-photon (2p2h7) contributions to the Z® self-energy in nuclear 
matter. The black dots represent Z^N —>■ yN elementary amplitudes and the dash lines denote strong interactions. To obtain the 
immaginary part, the intermedate states intersected by the dashed line have to be placed on the mass shell. 






explanations of the MiniBooNE excess of events bases on oscillations do not seem very plausible. Multinucleon 
mechanisms like those in Fig. 7 are being investigated. There are other possible explanations, like the possibility that 
a heavy neutrino is produced by weak [37] or electromagnetic [38] interactions in the detector, decaying radiatively 
afterwards. Such scenarios can be investigated in the forthcoming MicroBooNE experiment, capable of distinguishing 
photons from electrons. It should be recalled that any explanation of the MiniBooNE anomaly in terms of single 
photons, using the physics of the Standard Model or beyond it, should satisfy the NOMAD upper limit. 
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